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Abstract 

Winter oxygen depletion rates (WODR) (g O2°-m~?+d-') 
were determined for 13 lakes in central Alberta during the 
winter of 1982-83. Although oxygen decreased in all lakes 
for the first three and a half months after freeze-up, the 
decreases were nonlinear. The highest WODR were observed 
just after freeze-up. The nonlinear WODR were significantly 
correlated with estimates of lake productivity (total 
phosphorus (TP)), but were not significantly correlated with 
morphometry (eg. mean depth, P>0.20). 

When the WODR from the Alberta lakes were treated as 
linear, to enable a comparison with other studies dealing 
with WODR, a correlation was found between WODR and both 
morphometry and estimates of summer productivity. This 
relationship was significantly different from what previous 
investigators observed. When data from three other sets of 
temperate zone lakes were combined with data from this 
study, WODR were best predicted from a combination of mean 
depth (z in m) and mean summer TP (TPsu in mg-m-?) in the 
euphotic zone: 

WODRS= 04,0 ee 00247) TRsu ta. 0 \3arZ p= 0.90 
The above equation permits the prediction of WODR for a 
greater range of lake types than previous models. 

Models to predict WODR in lakes are based on oxygen 
profiles obtained from the deepest site in the lake. When 
the average dissolved oxygen concentration in Wizard Lake 


was calculated from oxygen profiles obtained from six sites, 
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the oxygen concentration was 47% higher than when the 
average concentration was calculated based on the main 
Sampling site. Thus, it appears that one-site sampling may 
not yield an accurate estimate of the winter oxygen content 


of lakes. 
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I. General Introduction 


"A skillful limnologist can probably learn more about the 
nature of a lake from a series of oxygen determinations than 
from any other kind of chemical data" 

CHutehinson, 91957, =p575) 

Dissolved oxygen (DO) in lakes is used by plants, 
animals, and bacteria for respiration. If the consumption of 
oxygen in a lake exceeds the input, a net oxygen deficit 
occurs. Most oxygen depletion studies concentrate on 
Stratified lakes during the summer. Thus most of the 
information regarding oxygen depletion in lakes is based on 
Summer oxygen data. To obtain accurate estimates of factors 
influencing oxygen depletion it is often necessary to bring 
tiogilake™ “intosthe Vaboratoryae'rnus avloteorethe 
information available on the mechanisms involved in oxygen 
depletion is from laboratory experiments using sediment 
eoresrand water from the laké. le is often diitiacule to 
extrapolate directly from observations of laboratory 
experiments to the situation in the field (Riley and Prepas 
1984). However, the general trends from the laboratory may 
be used to interpret the results from the field. 

In a stratified lake during the summer, DO is consumed 
throughout the water column. In the epilimnion, water 
fequlanly somes Winecontact with the airy thus DOvlevels ere 
usually close to 100% saturation. Water in the hypolimnion 


does not regularly come into contact with the atmosphere. 


Consequently, oxygen respired in the hypolimnion is often 
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not replenished and a deficit may develop (Thienemann 1928). 

Oxygen depletion in lakes is expressed in one of two 
ways: (1) volumetrically as mass of oxygen used per volume 
OFgwateragpernidaye:(qsma-ad=1)) or 12) -areallyoas imass sot 
oxygen used per unit surface area per day (g-m~?-d~'). 
Oxygen depletion rates were first used to define lake 
trophic status (Strom 1931; Hutchinson 1938). Recently 
Charlton (1980) and Cornett (1981) have shown that oxygen 
depletion is related to morphometry. Thus the earlier 
attempts by Strom (1931) and Hutchinson (1938) to classify 
lake trophic status based on the areal hypolimnetic oxygen 
deficit (AHOD) were incorrect. Oxygen depletion in lakes in 
winter 1S expressed areally since the primary site of oxygen 
consumption is at the sediment-water interface (Hayes and 
MacAuley 1959; Hargrave 1971) which is a function of lake 
area. 

An ice-covered lake in the winter is Similar to the 
hypolimnion of a Stratified lake. Both are isolated from the 
atmosphere and can develop oxygen deficits (Krumholtz and 
Cole 1959; Linsey 1981). With fall freeze-over, wind induced 
circulation ceases in a lake and atmospheric inputs of 
oxygen stop. Under ice-cover, DO is added during 
photosynthesis (Rhode 1955; Wright 1964; Jackson 1979) and 
by inflowing streams (Greenbank 1945; Pennak 1968). 
Concurrently, DO is lost through white ice formation, 
respiration and the oxidation of chemical compounds. When DO 


levels decline during winter, fish populations within the 


baadtel a jotrene en . o 
Desadacncs nepyxo 40" adam #4 yilnoit: mde br ; 
. meee ew Yiigasn’ is} sont ge ag - 


ae ES 
an ae 
A ee ‘Smee! qeb iq e276 soeiive gine 391 re xo 


eéeb oniiab os Heed Iesit sie. ester recat aie 
gidretag . (8021 ncaaisodgult 17 ee! morse ruled pangosd 
nepyko: eis avo avec! (TReT4 -saep409 bas (ORT t root yada: 


we 


poze) y. 
ms 7” 


a 


2 


0 


Ssifues odd aud? es Serags t6H q) batebeer as we casige®) . 
Wicseeéic oa beset) Hoenikawwietds CRT) moa qd esQastda” 
ABayNo otsenMtlegyd dagse end ad beaten edisde siagesd adel: _ 
fi wetel oc). neiseabgst aguas. 225954°95! S23W (COHAD sishish a 
genyie Jo atte yraming adj onal4 vifesis S8ecetgem BE —o - 
Gne 2ayer) eawdivasoi asite-wianties ands 26 


eiel it@ noi sondi 6 2: dospiw «98! -SverpssF 


S450 xeiiaje =i decciw saz. oi siel. Sesseoo-ei aa 
ans mat tietalee: sss dfod .sish beplisscte « Jo moa 
brid asionausd)) svitoiish nse7ed ga! camab 44 
bewinal. brew ,aevoresess) iis? 42th cAI GEr yeandsl seat 
72 2tuani Se teettg Sarre ene sdk Soi aoe Ger COE 

put o6 bends 2t of /e0oo-Ssei >45hnl cg 
bre (272) Acedcal s680). tigéae 3882) ehOGKI 4 
oe f8ART deames 5251 be jae 2m687 48. 
, AOLARNIGS del saidy decoy ape! 


oa cad atinscomis“okoags neta 


lake are subject to winterkill (Greenbank 1945; Scidmore 
1957; Halsey 1968+ Pennak 1968; Schneberger 1970; Casselman 
and Harvey 1975). Consequently, biologists are interested in 
predicting winter oxygen depletion rates (WODR) to determine 


which lakes may winterkill. 


A. Oxygen Sinks Under Ice-cover 

There are three major sinks for oxygen in lakes under 
1ce-coverz= (4) slushing, (2) achemical oxidation of reduced 
compounds, (3) biological consumption. These sinks are 
discussed in detail below. 

When black ice cracks, water seeps up over the ice 
flooding the snow and white ice forms (Adams 1976), this is 
termed slushing. As a result of this flooding, DO in the 
water which forms the white ice is removed from the lake. 
The total amount of oxygen removed from the lake during 
slushing events is proportional to the depth of white ice 
formed and the DO concentration in the flooding water. 

Chemical oxidation occurs when reduced compounds are 
Oxidized in theswater).columnsde.g. ,°2Fe.0,. +) 0295 2Fe20;. 
Reduced compounds (e.g., manganese, iron and sulfur) are 
released into the water overlying the sediments when DO 
levels are low (Mortimer 1971). These reduced compounds are 
transported via currents: and diffusion towards the lake 
surface and, if oxygen is encountered, they are oxidized. 


Brewer et al. (1977) found that the oxygen comsumption of 
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sediment cores was reduced by 91 to 99% by poisoning the 
sediments; the remaining 1 to 9% of the oxygen consumption 
was presumably due to chemical oxidation. Hargrave (1972) 
found that chemical oxidation accounted for only 20% of the 
oxygen consumption of sediment cores in the laboratory. 
Thus, as the amount of oxidizeable material increases, the 
consumption of oxygen by chemical oxidation increases. 
However, chemical oxygen depletion, in general, is not the 
main sink for oxygen in lakes. 

The major sink for DO is respiration by organisms 
(Hargrave 1969) including fish, invertebrates and microbes 
(egy fungi -and=bacterta). "Respiration by fish®accounts for a 
relatively small amount of the oxygen consumed ina lake. In 
Sharpe Bay (Jack Lake, Ontario), Linsey (1981) calculated 
EhaveeheaDOmresprred by fish (4.4 X40 mq-b.**d=*), was 
less than 0.2% of the total oxygen consumed. Benthic 
invertebrates also account for a relatively minor portion of 
total oxygen consumption in lakes. However, benthic 
invertebrates can greatly influence the oxygen uptake by 
bacteria (Edwards and Rolly 1965; Graneli 1979). By 
burrowing in the sediments and increasing the sediment 
surface area, benthic invertebrates increase oxygen 
consumption by bacteria. Benthic invertebrates also bring 
organic material and reduced substances to the sediment 
surface where they may be oxidized. Bacteria are the main 
consumers of DO in lakes (Edberg and Hofsten 1973; Brewer et 


al. 1977). A correlation was noted between bacteria 
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abundance and oxygen depletion (Zobell and Stodler 1940; 
hayes jand (“MacAuleyl959; eBrewer set 'alhor1977)2) Fungi play ia 
major role in oxygen consumption as they are responsible for 
the breakdown of large particles of organic matter (eg. dead 
macrophytes) into smaller particles suitable for bacterial 
consumption. Bacteria use oxygen for the decomposition of 
organic matter. AS more organic material is made available 
DOR Cecomposi tron py tthewnumber sof @baqteriasimereases and the 
oxygen demand by the bacteria inceases correspondingly 
(Zobell and Stodler 1940; Hayes and MacAuley 1959; Hayes and 
Anthony 1959). The site of highest bacterial numbers in 
lakes is at the sediment-water interface (Hayes and Anthony 
1959), thus the role of sediments is an important factor in 
oxygen depletion in lakes and will be discussed in the 


following section. 


B. Sites for Oxygen Depletion 

When considering oxygen depletion, a lake may be 
divided into two distinct compartments: (1)sediment-water 
interface and (2) open water. Oxygen depletion in these two 
compartments is termed sediment oxygen demand (SOD) and 
water oxygen demand (WOD), respectively. In deep lakes, 
Organic material falling out of the trophogenic zone is 
oxidized to a greater extent in the open water than the same 
material would be in a shallow lake, since the material has 


a greater distance to fall in the deep lake, hence a longer 
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residence. time, (Charlton, 1980;.Cornett, 1981)..1n laboratory 
experiments carried out on sediment cores and the overlying 
Water, =sEdberg- and Hofsten.(19/3.).found: WOD, accounted for.12% 
of the oxygen consumption when the water overlying the 
sediment cores was rich in oxidizable substrates; however, 
when water overlying the cores was from the outflow of a 
papermill, 40% of the total oxygen consumption was caused by 
WOD. AS the concentration of oxidizable material in the 
water column increases, the relative contribution of WOD to 
the total oxygen consumption in the lake increases. 

In lakes, the main site for oxygen depletion is at the 
sediment-water interface. Productivity and the mass of 
oxidizable material are relatively low in the water column 
during winter, therefore the contribution of WOD to the 
total oxygen consumption in the lake is small. In two 
Sub-arctic. Jakes, Chenard,(1980) found that SOD could 
account for between 50 and 100% of the total oxygen demand 
in the lakes during the winter. Hargrave (1971) found that 
oxygen consumption under ice-cover in Lake Esrom, Denmark, 
Gould be attributed almost entirely to oxygem uptake by the 
sediments. The nature of the sediments influences the rate 
of oxygen depletion and will be discussed in the following 


Section). 
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C. Factors Controlling Oxygen Depletion in Lakes 

Oxygen depletion rates in lakes are governed by two 
factors: ml )ethemnaturesand, amount=orethe: supstrates® too be 
oxidized (2) the supply of oxygen to the sites of oxygen 
consumption. In temperate zone lakes, most of the production 
of organic material takes place during the open-water season 
(Schindler and Nighswander 1970). Material produced in the 
open=water Season 1s partially: oxidized as. it» falls to the 
sediments where it is stored. Production in lakes rapidly 
declines when lakes freeze over. Thus, most of the material 
that will be oxidized during the winter is accumulated in 
the sediments from production that occurred during the 
ice-free season. Oxygen demand may increase as a lake 
undergoes eutrophication because of increased production 
(Brewer et al. 1977). However, Hargrave (1975) and Graneli 
(1978) found little change in oxygen uptake rate of sediment 
cores when fresh organic material was added. No correlation 
was found between the organic content of sediments and 
oxygen depletion rate (Edberg and Hofsten 1973). However, 
the nature and quality of the organic material in the 
sediments was not determined in these studies. This can be 
important since organic material can be composed of easily 
oxidizable material such as simple carbohydrates or material 
which is relatively difficult to breakdown such as lignins 
or humic complexes (Hargrave 1971). 

In laboratory experiments, oxygen depletion rates were 


higher at the beginning of the experiment than at the end 
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(Zobell and Stodler 1940). Zobell and Stodler (1940) 
attribute the reduction in the rate of oxygen consumption to 
a decrease in the availability of substrates suitable for 
oxidation. 

Rate of oxygen depletion is governed in part by the 
supply of oxygen to the sediments. Since sediments are the 
primary site for oxygen consumption, water in contact with 
the sediments quickly goes anoxic. However, as there is 
Circulation under ice-cover (Likens and Hasler 1962; Likens 
and Ragotzkie 1966) water is transported to the 
sediment-water interface, bringing DO to the sediments. 
Horizontal mixing is faster than vertical mixing in the open 
water. Vertical mixing in the open water is supplemented by 
the movement of water down the sides of the lake at the 
sediment surface. AS water approaches the sediment-water 
interface, it warms and becomes denser than the colder 
Surrounding water, and then slides down the sides of the 
lake displacing the bottom water upwards. This process 
results in a slow circulation of the lake water and brings 


oxygen to the sediments. 


D. Production of Oxygen Under Ice-Cover 

If no water enters the lake in winter, photosynthesis 
isetheronly source tof toxygen gaingin lakes. Lf snow cover on 
the lake is thin, enough light may enter the lake to permit 


a substantial amount of photosynthesis to take place. The 
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depth of snow needed to reduce light penetration to 
effectively halt photosynthesis would depend on the type and 
pattern of snow cover present on the lake and also the light 
regime of the lake (Adams and Roulet 1980). Barica et al. 
(1983) were unsuccessful in their attempt to increase 
photosynthesis under the ice by removing the snow cover. The 
presence of algae in the water column does not mean that 
there 1s oxygen production taking place in the water. 
Schindler and Nighswander (1970) found that the algal 
population under the ice in Clear Lake, Ontario, was 
relatively large but the algae were in a state of dormancy. 
Chlorophyll a (Chl a) levels in lakes, a measure of algal 
production, are relatively high just after freeze-up and 
decline throughout the winter (Greenbank 1945; Barica 1977; 
Jackson 1979). As a result, oxygen production caused by 
photosynthesis is negligible in most lakes throughout the 


the period of ice cover (Jackson 1979). 


E. Existing Models to Predict Winter Oxygen Depletion Rates 

Whole-lake WODR have been studied and modelled in two 
distinct lake types in the temperate zone lakes of North 
America. Deep oligotrophic lakes on the Precambrian shield 
have been studied by Welch et al. (1976) and the shallow 
eutrophic prairie pothole lakes have been studied by Barica 
and Mathias (1979). Both studies reported that WODR (g 


O,-m-?-d-') could best be predicted by models that used 
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morphometry (mean depth (z)) as the independent parameter. 
The relationship between z and WODR in the two studies was 
quite different: 


WODR OCOCR ORO 2 "ce (Welcheebrater 1o7Ge 


WODR = 0.14 + 0.062 z (Barica and Mathias 1979) 
Mathias and Barica (1980) suggest that this difference is a 
BUNGILON TOL productivityebut eheys-didinomrattempt. to 
incorporate an estimate of lake productivity into a model. 
Schindler (1971) related winter oxygen depletion to the 
amount Of allochthonous material received by the lakes in 
the Experimental Lakes Area (ELA). This relationship has 
never been confirmed in subsequent studies, probably because 
the vegetation cover in the lake basins in the other studies 
differs from lake to lake, whereas in the ELA the vegetation 
cover iS similar between lakes. Hence the source, type, and 
amount of allochthonous material entering the lakes is 
different for lakes in other regions (Mathias and Barica 
1930)" 

Jackson and Lasenby (1982) developed two models to 
predict oxygen profiles under 1ce-covered lakes. The two 
models are based on data collected for lakes on the 
Precambrian Shield and lakes in limestone basins in Ontario. 
The more productive limestone basin lakes had higher oxygen 
depletionerates, sndicating that oxygen depletion is) relaved 
to productivity. 

All of the individual studies on WODR have focused on 


groups of lakes with similar morphometries or 
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productivities. Thus, models based on these studies are only 
applicable to the lake type they are based on. To link 
together the previous studies, I gathered data on WODR in 13 
Albertan lakes, and examined the effect of morphometry and 
estimates of productivity (spring and summer total 
phosphorus, summer Chl a, and loss on ignition of the 
sediments) on WODR in these lakes. Furthermore, data from 
the literature and this study were combined to produce a new 
model to predict WODR over a broader range of lake types 


than the two existing models. 
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II. Modelling Winter Oxygen Depletion Rates in Temperate 


Zone Lakes 


A. Introduction 

Several researchers have investigated patterns in 
winter oxygen depletion rates (WODR) (Welch et al. 1976; 
Barica and Mathias 1979), factors controlling WODR (eg. 
Mathias and Barica 1980) and ways to reduce WODR (Wirth 
1970; Barica et al. 1983). Recently, Jackson and Lasenby 
(1982) constructed two models to predict oxygen profiles 
under ice in lakes located on the Precambrian Shield and in 
limestone basins in Ontario. 

Winter oxygen depletion rates have been studied in two 
extreme lake types in temperate regions of Canada: the 
relatively deep, oligotrophic lakes on the Precambrian 
Shield (Welch et al. 1976) and the shallow, eutrophic, 
prairie pothole lakes in southwestern Manitoba (Barica and 
Mathias 1979) (mean depths (z) ranged from 4-27 and 2-4 m 
and mean summer chlorophyll a (Chl a) concentrations in the 
Suphor tee zone. ranged £rom= 1-2vand=5=102emg= ms ~5 
respectively, for the two lake types). Although z was the 
best predictor of WODR in both lake types, the mathematical 
relationships were different (Welch et al. 1976; Barica and 
Mathias 1979). Thus, empirical models to predict WODR 
Perbaineto relacively=specific lake types, in different 
regions. Winter oxygen depletion rates are also affected by 


Snow cover (Barica et al. 1983) and the amount of organic 
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material produced in a lake (Mathias and Barica 1980). 
However, neither snow cover nor estimates of lake 
productivity have been incorporated into models to predict 
WODR. In contrast, models to predict summer areal 
hypolimnetic oxygen deficits (AHOD) do incorporate measures 
of both morphometry and productivity (Cornett and Rigler 
ino, 19803 Charltonyu980). 

Winterkill is a common phenomenon in lakes in central 
Alberta. Thus, considerable interest exists in factors 
controlling WODR in these lakes. However, Albertan lakes are 
different from other temperate zone lakes studied in Canada; 
they may be as productive as the prairie pothole lakes, but 
are much deeper (Prepas and Trew 1983). To test whether 
either existing WODR models are applicable in Albertan 
lakes, data were collected on 13 lakes. These data were also 
used to construct an empirical model to predict WODR for 
Albertan lakes, and were combined with other data in the 
Incterature to test the contribution of Lakes productivity to 
WODR. 

To estimate WODR, water samples are collected over the 
deepest part and treated as representative of the whole lake 
(Welch et al. 1976; Barica and Mathias 1979; Jackson and 
Lasenby 1982). This approach assumes that oxygen 
concentrations are homogeneous within horizontal strata; an 
hypothesis which has never been rigorously tested. To 
determine if one-site-sampling is adequate for estimating 


whole-lake oxygen mass, samples were collected from six to 
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seven sites in each of five lakes for comparison with 


resultsiei rom a@ Single ystation, 


B. The Study Lakes 

The “13 lakes sin this :study are wocated in central 
Alberta withinea 160ckm radius of ‘the city «of Edmonton.«The 
lakes range in size from 0.08 to 4.4 km? (see appendix B for 
bathymetric maps). All lakes are underlain by 10 to 40 m of 
Glacialstrll. Total sdissolvedtsolids in-the lakes ranges 
Prom fo0to) 220 ing*he Vand colour sanges Erom 2 s0025 mg:h 
Pt. The lakes were chosen to encompass large ranges in 
morphometry and productivity: z of the lakes ranges from 3 
womNOrn and trophic: Status vanreswirom oligotrophic sto 
hypereutrophic (mean summer euphotic zone Chl @ ranges from 
2 to 155 mg-m-*). Background data for the lakes are in 
Prepas and Trew (1983) and Prepas and Wisheu (1984). Two 
lakes, Amisk and Baptiste, have distinct north and south 
basins. Both basins of Amisk were sampled and only the deep 
south basin of Baptiste was sampled; all three basins were 
treated as separate lakes. Each of the remaining lakes was 
sampled at only one site. Seven of these lakes have simple 
Basins with a single deep spot, two lakes, Eden -and Sauer, 
each have two distinct deep spots and the remaining lake, 
Hubbles, has three deep spots. None of the study lakes have 
major inflows or outflows (Prepas 1983, unpublished), hence 
the main source of organic matter in the lakes is from 
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C. Materials and Methods 

The study lakes were permanently ice-covered by 5-11 
November, 1982, and remained frozen until the third week in 
April, 1983. Bach lake was visited seven to 13 times from 
Just prior to freeze-up until one month before break-up. On 
each visit, water samples were collected over the deepest 
part of the lake with an aluminum drop-sleeve water bottle. 
Water for dissolved oxygen (DO) analyses was collected every 
1 or 2 m from the top of the hydrostatic water level to the 
lake bottom or to where it became anoxic as evidenced by the 
smell of hydrogen sulfide. These samples were fixed 
immediately (Carpenter 1965) and, to avoid freezing, 
transported back to the laboratory in insulated boxes 
containing hot-water bottles. 

To determine the potential contribution of algae to the 
lake oxygen mass under ice-cover, water was also collected 
for Chl @ analyses. For the first one and a half months, 
samples were collected at i-m intervals from 1 to 4 m and 
for the remainder of the winter, samples were collected to a 
depeneot Seto 24ym and@pocledsintorstrata 2 to Gi medeep. 
Water samples for Chl a analyses were placed in 2-L amber 
Nalgene bottles and transported to the laboratory. To 
estimate spring and summer productivity in the lakes, total 
phosphorus (@P) and$¢hiea levels»ingthe euphotre zone for 
ten of the lakes are from Prepas and Vickery (1984). For two 
lakes, Nakamun and Halfmoon, TP and Chl a@ values are from 


Riley (1983). I sampled the 13th lake, Hubbles, five times 
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fromtMay*i4%to Auge 26,;1)1982: <In«this study, water samples 
for TP and Chl @ analyses were collected with Tygon tubing 
at three stations and pooled. These samples were collected 
Feomethe euphotic zone, idefinedeas thesdepth=to  1%ssurface 
Peradianece, of the vphotosyntheticalliy active: radiation, as 
measured with a Lambda L1-185 light meter equipped with an 
underwater sensor, or 2 times Secchi disc depth (SD). 

To determine whether organic content of the sediments 
1S an estimate of lake productivity which can be used to 
predict WODR, sediment samples were collected from 12 of the 
study lakes during the summers of 1982 and 1983. The number 
of samples collected (4-8) on each lake increased in 
proportion to lake size. These samples were collected 
randomly from each lake with a four-barrel corer (Hamilton 
et al. 1970). The top 5 cm was removed from each core and 
pooled with other cores from that site. The samples were 
broughteback tothe baboratory “and frozensuntrl they could 
bewanalyzed for Loss son ignition (LO). 

Detailed morphometric data were already available for 
five of the study lakes (Prepas and Trew 1983, unpublished). 
The remaining eight lakes were sounded with a Furuno model 
Fe-400 depth sounder. Morphometric data were determined from 
thembathymetric maps with ‘a Tallos digitizer connected toa 
Hewlett-Packard 9825B desk-top computer supplied by Alberta 
Energy and Natural Resources. 

Dissolved oxygen concentration was determined on 


duplicate water samples based on Carpenter's (1965) 
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modification of the Winkler technique. Most samples were 
titrated within 12 h of collection and all were titrated 
within 24 h. To determine the accuracy of the method for 
collecting DO samples, I compared the average DO values from 
duplicate water samples. The average difference between 
independent oxygen samples collected at the same site and 
depths was0207eme<k *¥and. the range was 080d tor0.322omg-L-". 
Percent oxygen Saturation of the lakes was determined based 
on oxygen saturation tables in Cole (1983). Duplicate 
samples for Chl @ analyses were prepared within 12h of 
eollection by filtering’ 0705° to 2 L of lake water, edepending 
onmeec lls density, through glass. tiber ta lters) ehiorophyl1 +a 
was determined with the ethanol extraction technique of 
Bergmann and Peters (1980). Total phosphorus was determined 
on duplicate 50 mL samples of lake water with the persulfate 
method of Prepas and Rigler (1982). Loss on ignition (LOI) 
was determined by the weight loss of dried sediment samples 
baked at) 550 C for 24)h. The’ LOI was not corrected for 
dehydration of the samples (Hargrave 1975). Loss on ignition 
estimates for all samples on each lake were combined to 
yield an average whole-lake estimate of LOI. 

For each winter visit, oxygen mass was determined for 
each stratum from the product of stratum volume and oxygen 
concentration; these data were summed to yield whole-lake 
oxygen mass. The volume of the top i-m was corrected for ice 
thickness. The whole-lake oxygen mass was divided by the 


under-ice surface area of the lake to yield the areal oxygen 
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mass (g O2*m~?) (Welch 1974). To calculate WODR (g 
O2°-m-*+d~'), areal oxygen mass was regressed against days 
past freeze-up; the slope of this regression line is the 
WODR. 

To develop empirical models for the lakes in central 
Alberta, WODR were regressed against estimates of spring and 
Summer productivity (TP; /Chlya, and. LOI)Vand mean and 
maximum depth (z and Zmax, respectively). Volumetric 
(mg-m-*) and areal (mg-m-?) expressions were used for TP and 
Chl a; the areal expressions were calculated as the product 
of TP and Chl a@ concentrations in the euphotic zone and the 
depth of the euphotic zone. To construct a general model to 
predict WODR, data (j.e.:, WODR, z, TP,.and,Chl a) from other 
temperate zone lakes (Schindler 1971; Dillon and Rigler 
1974: Reid et al. 1975: Welch et al. 1976; Barica et al. 
1978; Barica and Mathias 1979) were combined with data from 
this study and analysed by regression analyses. Spring was 
defined as 27 May-14 June for Alberta data collected in 
1982. Summer was defined as 15 June-15 Sept. for Alberta 
data collected in 1982 (after Prepas and Vickery 1984), and 
end of May until end of August for unpublished data 
collected on Ontario lakes. Spring and summer data from 
published sources were used as published. 

The regression analyses were performed with the BMDP 
Statistical package program PIR (Dixon 1981) on the Amdahl 
model 5860 computer at the University of Alberta. Other 


Statistical analyses were from Snedecor and Cochran (1980) 
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Or Prepas (1984) "All Walues in this report aré expressed as 
means + standard error, except where noted. 

To evaluate the differences between predictions of WODR 
with summer areal hypolimnetic oxygen depletion (AHOD) 
models and WODR models, I compared the residual mean squares 
(RMS) for each of four AHOD models with the RMS for the WODR 
model(s). The RMS for each model was calculated as follows: 

RMS -2 (pred,-obs,)*/df 
where pred; is the predicted WODR (from the AHOD model) and 
obs; 1s the observed WODR for each lake (i), and df is the 
degrees of freedom. For each model, df was defined as number 
of data points minus one minus the number of independent 
variables in the model. Since multiple comparisons were made 
in this analysis, the chance of committing a type I error 
was greater than the value given in a standard F-table; 
therefore, a probability level of 0.02 was used to 


distinguish significance. 


D. Results 


Fall Turnover 

Just prior to freeze-up, DO levels were surprisingly 
low in 12 of the 13 study lakes; the average oxygen 
Saturation level was calculated to be 69% at the time of 
freeze-up, with values as low as 51%. Only one lake, Hasse, 


was almost fully saturated (94%) (Table 1). These findings 


at eee . . 
GOW To notisibesG g2ev73d 
_ a 


(aGiA) agi sasqen nap sno Si genttt Loar 


assups 


ine nso beiiStzes ait feskquo> t ,aiesbow iow bie Blebor 4 


“ 
tessa brebal 

$60 
' i 


ea 
i a 
7 


Sal edn va tis ets 


io 


ty does tod (2 


ie 7 ook ee 
MON ed? tot SMA sda Gsiw 2acssem GGHA: 7407 
o 
seuclin® te ibesotivisles ef. [Sbam tocs ror ake om . (adtsbon 
: Lal ~ « —_ 7 
$842 ( eto- ,663q).— = =a 
7 
GG | labam doks. 40% ter) BION Helniheryg ads ecg erg sitode 
- ’ ‘ . q ’ 
sai at = DAY _f 7) ai4us° aes oi ee ie Hat 2 wT a¢ ‘ade 
s 7% tabo: - fi 4% 3a fs 29h 
newt 26 Healt mw 26 . stom dose . is — 
jaaboacsbar 26 Sadun #4 aveja-eho eure 2inicg east 
‘ 
shan grew efori tatters sigirivuem edt ce smo ens AL eolde 1s 
ne 
1975 T Sv? 6 pres? cine: afi cir Oye eile atay net 
i = i 
raider S262 ne84 it \ oft cats t8efnewg Bae 
ot bse 2e- 0 %o Levef era ciaedosg ,sootsitea2 
wires tote dabupnd za vi6 
<eremeen 
‘lenieiaqiusa ataw efevel O8% (Qb-eree33 o¢ TTI, gee 


WSeveo seez7eVAa ois 


7 2G. ont: ada ds ea: ad o 
_ 
i ve print? “eae 


ha 


Sita i 


F 


7 
al _ 


7 
re 


Saal om wes, giao. we 26 We br @q.a ae 
ati ou ee aes axusee gt 


‘eeeder yeuse 61 


gad poner 


4 nee eae peri a1 
os 


ee 


ai = 


ie 7 va 


ae 


are contrary to observations in other areas; e.g. in the 
pothole lakes of Manitoba, oxygen saturation levels are 
80-100% just prior to freeze-up (Barica and Mathias 1979). 
Jackson and Lasenby (1982) reported that the maximum 
volume-weighted oxygen content of the Ontario lakes that 
they studied, was achieved approximately 15 days after 
rreeze-upy thiswise notmthéucases inmthenAlbertangdakes.oin 
only four of the 13 lakes, the saturation level increased 
between the time of freeze-up and the first sampling date, 
possibly due to relatively high oxygen consumption relative 
to the inputs of oxygen via freezeout and photosynthesis. 
Not only were DO levels low in the study lakes in the fall, 
but mixing was incomplete in two of them. A layer of anoxic 
water up to 5-m high remained over at least 18% of the 
sediments in the deepest part of Twin Lake. In Hubbles Lake, 
up to 9-m of water, overlying 15% of the sediments, was 
anoxic. Since oxygen consumption in lakes is mainly caused 
by bacteria at the sediment-water interface (Brewer et al. 
1977), incomplete circulation could reduce oxygen depletion 
rates in the water column. Consequently, for the subsequent 
analyses, WODR for Twin and Hubbles lakes were not pooled 
with those from the remaining 11 study lakes where DO levels 
were above zero over the entire sediment area prior to 


freeze-up. 
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Lake Cover and Chlorophyll a 

During the visits to the lakes while they were 
1Ge-covered, @icetthicknessvat the sampling sites averaged 
44.1 + 1.35 cm and average snow depth was 8.2 + 0.94 em. 
White ice was noted only on Twin Lake; it was 2.5 cm deep. 
The total snowfall at Edmonton International Airport during 
the winter of 1982-83 was 0.92 m, which was lower (t-test, 
P<0.001) than the previous 30-yr average (1.25 m) 
(Environment Canada). Consequently, the amount of light 
entering the lakes during the winter of 1983 was greater 
than during an average winter. 

PnP i 2n0t (the tstudy glakesm Chi a concentravionsin the 
top 4 m decreased with the onset of ice-cover and remained 
low until late February or early March, as illustrated with 
data from Lessard Lake (Fig. 1). In these 12 lakes, average 
Chl @ concentration in the top 4 m between freeze-up and 
Marchoeco, 1983) (was 2 ..98t50755 mg=meewhile aitersMarch 691¢ 
rose to 4.6 + 1.42 mg-m'°*. For the same lakes, average snow 
depth during these two periods was 7.3 + 0.89 cm and 10.3 ¢ 
1:15 cm, respectively. In the 13th lake, Halfmoon, a bloom 
Stmwalgae (32 mg-m**GChi a) was recordedvat od depthwor 1 ymeon 
December 14 but decreased to 3 mg-m~* Chl a@ 10 days later. 
THES algae bloom could have intreduced oxygen under the ice 
but the DO level at 1 m decreased from 6.6 mg°L™* on 
December 148to0 496 "mg-L** on December 2¢.8Snow>cover von 
Halfmoon Lake for these same two dates was 4 cm and 5 cm, 


respectively. In March, 1983, DO and Chl @ levels increased 
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in the top few meters in 10 of the lakes as illustrated in 
Pigs. 2uvith datasfirom Sater Lake. InaSauerubake,eChlia 
merveased five folds (fromil3.3eto 15+Semg-m 2)einuthe tepe4 m 
whale DO increased 18 fold’ (frem-0506 to 1¢2emg-Bet) einsthe 
Same stratum between March 2 and March 25, 1983. Since DO 
and Chl a@ levels began systematic increases under ice only 
after early March,1983, photosynthesis probably contributed 
an insignificant amount to oxygen levels between freeze-up 


and early March, 1983. 


Dissolved Oxygen 

For the first 112 days past freeze-up, the mass of DO 
(g-m”-?) decreased systematically in all lakes. Beginning in 
early March, changes in the mass of DO followed two distinct 
patterns: in three lakes DO increased while in the remaining 
10 lakes it continued to decrease, but at a lower rate than 
previously. Since DO levels systematically decreased in all 
lakes during the first 112 days past freeze-up and during 
the same period there was no evidence that a significant 
amount of oxygen was added to the lakes via photosynthesis, 
WODR were calculated for this 112-day period. Winter oxygen 
depletion rates, although systematic, were not linear during 
the 112-day period; rather WODR were more rapid at the 
beginning than at the end of this period. To illustrate this 
pattermp, therdataston theoll2-day period were subdivided 
into two periods: (1) the first three sampling dates and (2) 


the remaining sampling dates, and WODR were calculated for 
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both periods. Since a minimum of three points are needed to 
construct a regression line, subdivided WODR were only 
calculated for the eight lakes which had a minimum of six 
Sampling dates during the 112-day period. For these eight 
lakes, the WODR for the first three sampling dates were 
consistantly higher than rates based on the remaining dates 
(WIteoxon Signéd-rank test, P<0.05) (Table w2)s These results 
were Surprising; WODR are traditionally considered as linear 
until average whole-lake DO levels fall below 3 mg-L-~' 
(Mathias and Barica 1980). My data do not support this 
hypothesis; e.g., in two lakes, Lessard and Eden, average 
whole lake DO levels remained above 3.6 and 4.0 mg-L"', 
respectively, during the 112-day period when WODR were 
calculated. Yet in these same lakes, WODR just after 
freeze-up were 54 and 68% higher, respectively, than later 
thether 12-day period (Table 2) eyThus; sini’ thts study, WODR 
were nonlinear regardless of DO levels. Since changes in the 
mass of DO were nonlinear in the study lakes, a suitable 
transformation was sought for these data. The best 
transformation for the mass of DO (X) vs. time was a power 
aineeion: XK “=]|-X%"9 *?. For 12 of the 3137 lakes, the 
correlation between the transformed DO data and days past 
freeze-up was better than or equal to the correlation based 
on the untransformed data (Table 2). These data are the 
first WODR treated as a nonlinear process. Thus these data 
cannot be compared with any existing WODR models. I examined 


the relationship between two morphometric parameters (z and 


wide sees +03: ere amie pi onan! ae 


f& 


: - wo 
\ anew dasabs pralemee asa03 32 ys! etd 203 AGGw 


f i. . aad f - : iy a 
espcveve. (s°bs O76 bY sada svt. y 2 
mal 
iden 0.2 ie 2,0 evoes Ssqvenai 2ievel VA oe 
: . 
ead pao nanw r eh 7 5 Sis. « art i) = yiewi 
Seite, taur Het ~eetsl sities! seen: ; ge7. 
> = = < fr | “Tie ao ¢ =36%8 
= od ot ee io ap | (=< ; 2 = int : 


fe 


e423 nl @egtatts s5nit jeiavel oT tc Saale; ape? cies 
el@gttca eeroded YOulp@-o01 lip, sHgatd oe - 

5 ve: 

faeiteh? Gaiso wesc: 253 TheUse ae nei 
" 


jewel a Seu eeit sax call getty box on c 

J 

7 J iy : Pi. 7 oY P 
an esist £7 ae) Shinai ee = 


a. t2pn Sveb ce Breb Go -Gense {2s 1t2) acs ae 


_ . Beaed sal i lek Bo oF Lees, Yo. 568s a 


wild e718 haps AenT 1! tS ei any a 


vail sata ters a 7 
a ; a 4 m 


Zo 


zmax), four parameters used as estimates of annual 
production of organic matter (LOI, spring TP, summer TP and 
summer Chl @ levels) and WODR based on the transformed data 
(WODR') (Table 3). Winter oxygen depletion rates based on 
transformed data (WODR') were positively correlated with all 
variables except LOI. The significant negative relationship 
between LOI and WODR' is hard to explain. Surprisingly, the 
only positive significant relationships were between 
estimates of open water productivity (TP and Chl a) and 
WODR' (Table 3); the best predictors of WODR' were summer TP 
andsChi av(TPsu and Chl 7a) ™intmo=m- 2: 


WODR" = 0700554 0.00047 TPst (ap) 


Ml 


WODR' OF062> #-0300052 Chia (2) 

These results are contrary to investigations on two other 
groups of lakes (Welch et al. 1976; Barica and Mathias 1979) 
where DO data were not transformed prior to calculating 
WODR. In the other investigations, morphometric parameters 
were much better predictors of WODR than estimates of summer 
productivity, whereas summer productivity was the best 
estimator of WODR with the transformed WODR data from this 
study. 

Although data from this study suggest that models to 
predict WODR are best based on transformed DO data, 
untransformed data must be used to compare WODR in Albertan 
lakes with temperate zone lakes in other regions. Thus the 


subsequent discussion is based on WODR calculated with the 


Ofiginal data. 
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In the 13 Albertan lakes, WODR covered twice the range, 
Ovi 2omtoR0'848 GiGs =m isd a fas thateobsenveaeim theatwo 
previouswsStudies on WODR jin Canadian lakes (Bigs 3)meThe 
WODR for the Albertan lakes were tested against existing 
empirical models. The model for lakes on the Precambrian 
Shield which are similar in depth but less productive (Welch 
et al. 1976) consistently underestimated WODR in the 
Albertan lakes (Wilcoxon signed-rank test, P<0.01). 
Conversely, the model for the generally shallower prairie 
pothole lakes (Barica and Mathias 1979) consistently 
overeStimated WODR in the Albertan lakes (Wilcoxon 
srqned rank ‘test, (P<0505) (Tablet4, Fig.63)% 

To compare WODR in Albertan lakes with the lakes on the 
Precambrian shield and in southern Manitoba, WODR in the 
Albertan lakes were regressed against parameters indicative 
of lake morphometry and productivity (Table 3). Based on 
untransformed data, patterns in the Albertan lakes were 
comparable with the studies in other regions (Welch et al. 
1976; Barica and Mathias 1979); morphometry was a stronger 
predictor of WODR than estimates of productivity. The best 
predictor of WODR was mean depth (z in m): 

WODRE=—0e20 +00. 0362 (3) 
Unfortunately, too few Albertan lakes were sampled to 
construct multiple regressions on these data alone. 
Productivity, as estimated by spring and summer TP, was 
significantly correlated with WODR only when TP was 


expressed on an areal basis. Thus all further calculations 
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are based on areal expressions of productivity estimates. 

A serious limitation with the existing models to 
predict WODR is that they work only in the region where they 
were developed. These published models to predict WODR are 
based only on morphometry (Welch et al. 1976; Barica and 
Mathias 1979). However, there are numerous suggestions that 
oxygen depletion rates are a function of both z and lake 
productivity (Tablie 3p ithis study; Mathiasewand Barica 1980; 
Charlton 1980; Cornett and Rigler 1980). Previously, data 
were insufficient to test whether a model incorporating both 
morphometry and estimates of productivity: (1) could be 
developed to predict WODR and (2) would be generally 
applicable to north temperate zone lakes. To develop a model 
to predict WODR over a wide range of lake types, estimates 
of WODR, open-water productivity, and morphometry from 48 
lakes in Alberta, central and northwestern Ontario, and 
Manitoba (Table 5) were combined. These lakes covered a 
broad range in terms "of mean depth “(1.5<'2 322-7 m), “and 
estimates of ‘summer productivity nes Chima <184 mg:mr~ and 
PiaeTP 4222 mg*m *) inthe euphotic zone. Untortunately a1P 
data were not collected for the Manitoban lakes and TP 
estimates were available for only six lakes in northwestern 
Ontario. In the combined data set, WODR were strongly 
correlated with estimates of lake productivity and 
morphometry (Table 6). However, contrary to the results 
reported for regional data sets (excluding the transformed 


data set for the Albertan lakes), estimates of productivity 
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were Stronger single predictors of WODR than morphometry. 
Thus morphometric variables are better predictors of WODR 
within regions of similar productivity, but between regions 
of uiteringiproductivaty, productivityarser hes strongest 
predictor of WODR. Next, I combined z with the estimates of 
productivity to examine whether a multiple regression would 
improve on the linear models to predict WODR. In all three 
cases, including estimates of both morphometry and 
open-water productivity improved the WODR models (Table 6). 
In each multiple regression, the estimate of open-water 
productivity explained more of the variation in WODR than 
morphometry. The WODR model based on summer TP (mg-m~?) and 
z had the highest correlation: 

WODR* =) 02.1 010 450).00247-TPst) +02 013442 (4) 
However, TP data were available for fewer lakes than Chl a. 
When regressions were run on the lakes where summer TP data 
Were avallables) buEyinceorporatang WODR, Chlea, andi2zvin the 
regression, the correlation was highly significant (r = 
0.78, P<0.0001) but lower than when summer TP was used. 
Thus, summer TP is a better predictor of WODR than Chl 4a. 
However, models based on z and some estimate of 
Peoddetivaty sy Teeorechina, Can be usSed!l toppredicte WODR over 


aewide range of lakes in north temperate regions. 


One-Site-Sampling 
In the five lakes where DO profiles were taken at more 


than one site, considerable horizontal variations Ln DO 
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levels*were -found /*There “was “noeconsistent pattern in=the DO 
profiles at the various stations; depth at site sampled did 
noteingeluence WO "levels. For “anstanceysone22 January, 1983, 
in Lessard Lake at a depth of 1 m, DO concentration 
piereased from ™7.4 mg:b9" at the site of maximum depth to 
9.6 mg-L~' at the site where the depth was only 1 m; whereas 
one o-January, 1983; “in=Halfimoom Lake ata depth ot "lim, DO 
concentration decreased from 2.1 mg-L~' at the site of 
maximum depth to 0.9 mg-:L~'’ at the site where the depth was 
only 1m. To illustrate the possible error when only one 
Site 1S sampled and the whole-lake oxygen mass is based on 
that site, oxygen profiles collected at six sites in Wizard 
Bake“on 7 March, 1983, (Table 7) were used: To calculate DO 
mass for the whole lake based on all the sampling sites, DO 
mass was calculated for each site and summed together. To 
calculate the mass of oxygen at each site, the volume of 
each stratum was divided by the number of samples taken at 
that stratum depth. The DO concentration of each sample from 
each site, was multiplied by the corresponding volume 
associated with the sample, and all were summed together to 
provide the whole-lake oxygen mass. Whole-lake oxygen mass 
was divided by the volume of the lake to yield average 
oxygen concentration. The average DO value was 47% higher 
when data from all stations were used rather than data from 
Siees ibe wor maximumedepten (252 “and ).5eng-Ly 7, 


respectively). 
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To further illustrate the possible errors which could 
result when whole-lake oxygen mass is calculated from data 
collected at only one site, three separate estimates of 
oxygen mass (g-m~?) were made for the upper 8 m of Wizard 
Lake. At each station that was a minimum of 8 m deep, the DO 
concentration at each depth was multiplied by the whole-lake 
stratum volume corresponding to the depth of the sample. The 
results were summed to yield oxygen mass in the upper 8 m 
and divided by the under-ice surface area. Oxygen mass 
caleulated in thistmanner manged fromsécittorl6/8 “q-m-+2 "the 
Maan estatron had 18.35 q*m=*. sThus, “it appears that sampling 
at more than one site would improve estimates of whole-lake 
oxygen mass over one-Ssite-sampling for Albertan lakes. 
However, in all sampling programmes, a trade-off must be 
made between the number of lakes sampled during the winter 


and the number of sites sampled on a lake at each visit. 


E. Discussion 

Previous authors (Welch et al. 1976; Barica and Mathias 
1979: Mathias and Barica 1980) have reported that WODR are 
linear throughout the winter as long as average lake DO 
remains above anywhere from 1 to 3° mg*b™4. Data from this 
study do NOt Support ithe generality of (this observat ton sOne 
possible explanation for the non-linear trends in my data is 
that WODR for the Albertan lakes are based on a longer time 
interval than previous studies. The WODR for the prairie 


pothole lakes were based on approximately 90 d (Barica and 
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Mathias 1979) as were the Ontario lakes (Welch et al. 1976). 
The nonlinear trends in WODR in this study can be explained 
if some organic material in the sediments is more easily 
oxidized and is used up first (Edberg 1976; Jackson and 
Lasenby 1982). The rate of oxygen supply to the sediments 
may also influence WODR. Since oxygen depletion rates are 
greatest at the sediment-water interface, whole-lake oxygen 
depletion is, in part, governed by the rate of oxygen supply 
to the sediments. The rate of oxygen supply to the sediment 
is a function of the distance to the sediments from water 
containing oxygen and the concentration of the water 
containing oxygen: 

oxygen flux rate = K*(oxygen concentration gradient) 
where K is the eddy diffusion constant (Mathias and Barica 
1980). All of the Albertan lakes that mixed to the bottom at 
fall overturn developed clinograde oxygen profiles during 
the winter as shown by data from Baptiste South (Fig. 4). As 
the winter progressed, the DO in the lake decreased, the 
distance from the sediments to water containing DO 
increased, thus the oxygen flux rate could have also 
decreased. Methods to accurately measure the movement of 
oxygen from open water to the sediment-water interface need 
to be developed to assess the effect(s) of changing oxygen 
supply rates on WODR. 

To test whether incomplete mixing at fall overturn 
would depress WODR, the data on Hubbles and Twin lakes were 


examined. Total phosphorus in Twin Lake (69 mg-m™*) is 


ceri 


* Garer 4 a) api al) 


sigue ‘ee na: re 
yiiess. sion ahem each 7 aed nt eee 


a 


Hrs 


a“ lout bee 


ASP {AC 


\ lees 


*aemiboe Strz 


nbeaosect 


° 
SAGs -£< , 


Agpyas 


— 


7 - 
Ls Riiars ee 2 ; 
bat piuagor : f= a 


a - 


wit 7: vleota: fspyeo +o 


‘opislgep aabeeo Sante |. 4408 eyi(eird 2 


(aay  shsgzaeiit Ts2a” Ireamroee aaa 
so @2b7. Ht Od MeAgsvoOn:. Seq. peer 
x mi 
- ie ? - sAt = z “ « ya 
4 yigve, FRBYAU ay Sts rf onekbitea 
Stnemibes set o2~-22ne72 9 i425 ae 


eH4 2¢no iT En nese) 34 (O74, Nepyee, pith ee 


neil 


(snetbat aciters Js es + 5 fi nate - 
galrek tas 26 3M oe , Wiss ® yibs sds er) a 
i ae 
& MO9s0q si- of fii. 2en? Ss : +. 254, 2 2a, 
nn i eentée font MSC ¥eO aoetTooni 2 Bi Sas ee 
Phu See te r : 
' 
sal 208: (pie) Aaude ose feted’ nor? $388 qd’ nwone ae eee 
a> Sane 751 + i : ior 7 2 oe By ig t 34a ps — * 
OG Priiniesiqd> We3sy' Os ike te Nn ‘pod errr 
é 
a a . ; i ‘ 2 
oelé evad Giliod “ses xia Rebirs af?’ eifi?'s 
in. shomeyoniat} \escges..7fsPa9oI2e OF abe digecye 
i = LL el ee Pog ar eo 7 z 
> : v4 = <4 +s | Paes +» iar 
* ree 9 - THs TS se0-inen2 96 sha Ie 3 1534200 ayo 


natve6) od ipdans So, la) Saetie edt. angus 0) 


aciaissgnto sa 25> ont bt ant Ba 


- aoe sonst 


a 


af m + pit ea} $489) nh A 
ay 


ett ad nie te 


e 


nt rant ss 
nwt bap olagin, no. - 


aD 


twit 
- 7 - 
a7 


- 
a 


7 


36 


outside of the range of TP values used to construct hep 
(97 to 223 mg-m-?), thus no comparison should be made 
between the predicted and observed WODR' for Twin Lake. The 
observed transformed WODR for Hubbles Lake (0.0846) was 
inside the 95% confidence interval predicted from Eq. 1 
(0.0576 + 0.0472), however the observed value was 41% higher 
than the predicted value. A possible explanation for why the 
observed WODR' is higher, rather than lower, than the 
predicted rate is that Hubbles Lake was sampled at only one 
location and there are three distinct deep spots, thus the 
error in the measured WODR may be greater than in the other 
Albertan lakes. Both Twin and Hubbles lakes have summer Chl 
@ levels (6 and 21 mg-m-*) below the range of Chl a values 
Heed to generate™ig. 2 (27°%o° 147 mo-m- 2) ,-thus ac is? not 
possible to accurately predict WODR from Eq. 2 for these two 
lakes. Both Twin and Hubbles lakes have summer TP and mean 
depth values within the range of data used to generate Eq. 
4-)the predicted WODR, calculated from Bq. 4, for these two 
lakes (+ 95% confidence interval) were: Hubbles, 0.311 + 
Omiocseand Twine. 0se7 0g Os °m> == de Oey OP W0Z ei hepopserved 
WODR for the lakes were: Hubbles, 0.475 and Twin, 0.193 g 
O.°-m-?-d-‘. As noted earlier, the observed WODR for Hubbles 
Lakes is probably incorrect due to inadequate sampling. 
Thus, more data are needed to determine the effect(s) of 
incomplete fall overturn on WODR. 

To obtain better estimates of the oxygen content in 


lakes, data from this study suggest samples should be 
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collected at several stations. Further research is needed to 
determine the extent, magnitude and cause(s) of spatial 
heterogeneity of DO profiles in lakes under ice-cover. Some 
lakely causes are: (1) groundwater entering the lake, (2) 
differential oxygen depletion at various sites in the lake 
(Linsey 1981) and, (3) localized inputs of oxygen to the 
lake. 

Jackson and Lasenby (1982) constructed two models that 
predict oxygen profiles in lakes under ice-cover: one model 
for lakes situated on the Precambrian Shield and another for 
lakes situated in limestone basins. These models are based 
Onethe oxygen concentration an theilake Seto 15 .dsatter 
freeze-up. The two models were tested against the Albertan 
lake data to determine if oxygen profiles could be 
predicted. The model developed for lakes located on the 
Precambrian Shield consistently overestimated the final DO 
levels in the Albertan lakes by an average of 2.5 mg°L™'; 
conversely, the model for lakes in limestone basins 
consistently underestimated the final oxygen levels by an 
average of almost 2.0 mg-L~'. An underlying assumption of 
the models is that lakes in limestone basins are more 
productive than those on the Precambrian Shield; hence, DO 
is depleted faster in the limestone basin lakes. Jackson and 
Lasenby!s model for lakes on the Precambrian) Shield 1sabaced 
on only five lakes where the maximum depth ranges from eight 
to 33 m and average summer Chl a@ is approximately 2 mg-m°*. 


The model developed for lakes in limestone basins is based 
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on only three lakes where the maximum depth ranges from 16 
roe To me TSUuMMe RO rOduchIVvaty estimates are mGmraval babae for 
these lakes, however they are probably more productive than 
Ene (shield lakes but less productive than the Albertan 
lakes. Since the Albertan lakes are more productive than 
Ontario lakes, both “on the Precambrian Shield and in 
limestone basins (mean summer Chl a concentration (mg-m~* in 
19 Ontario lakes was 1.4 (Dillon and Rigler 1974) and in 27 
Albertan lakes was 26.5 (Prepas and Trew 1983)) it was 
expected that both models would predict higher final oxygen 
levels than observed in the Albertan lakes. Thus, the 
results were unexpected. However, it is possible that with a 
larger data set, incorporating an estimate of productivity 
into models which predict oxygen profiles under ice-cover 
will improve their generality. 

If the same processes are influencing WODR and summer 
areal hypolimnetic oxygen deficits (AHOD), AHOD models may 
predict WODR. Four AHOD models were chosen and tested with 
WODR data from the Albertan lakes. The four AHOD models 
chosen require inputs which were available for the Albertan 
Vakes estimates Of -open-water sproductivuty (SD, summer tP 
and Chl a@ levels), morphometry (z), and water temperature 
(Gaserby 1975; Charlton 1980; Counert tandeRigler 91980 )iihe 
assess the accuracy with which the AHOD models predicted 
WODR, I compared the residual mean squares (RMS) from the 
predicted WODR and the observed WODR for the 11 Albertan 


lakes (Table 4) with the RMS from Eq. 3 (RMS of equation 3 
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is 0.0078). Although the RMS for the AHOD models were four 
to 50 times greater than for Eq. 3, only Lasenby's model 
explained statistically less variation than the WODR model 
on=the- Albertan” lakes* (F-test; two tailed, P<0.005).° The 
AHOD models by Cornett and Rigler, and Charlton each 
explained approximately the same amount of variation of WODR 
in the 11 Albertan lakes, and although the RMS's for the 
AHOD models were four times greater than the RMS from Eq. 3, 
the differences were not statistically significant (F-test, 
evo falled,” 0% 05=P<05200). 

To determine if the AHOD models could predict WODR over 
a broad range of lake types, I compared the RMS from Eq. 4 
(RMS = 0.0065) to the RMS for the two AHOD models of Cornett 
and Rigler, based on the same lakes used to construct Eq. 4. 
The RMS for model #3, Table 4, based on mean summer TP and 
mean depth, was 0.0124 and the RMS for model #4, Table 4, 
based on mean summer SD and mean depth, was 0.0137. Although 
the RMS for the two AHOD models are approximately twice as 
great as the RMS for Eq. 4, the differences are not 
Statistically#signitteant (0.02=P<0510)% Thus, at appears 
that WODR can be accurately predicted from equations that 
incorporate estimates of both productivity and morphometry. 

The most important parameter in the models to predict 
WODR and AHOD over a broad range of lake types are estimates 
of open-water productivity (summer TP or Chl a) (This study, 
Table 6; Cornett and Rigler 1980; Charlton 1980). Both WODR 


and AHOD are influenced by morphometry but to a lesser 
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degree than productivity. Summer hypolimnetic oxygen 
depletion is also influenced by temperature (Charlton 1980; 
Cornett 1981) and WODR could be as well, although sufficient 
data do not exist to test this hypothesis. Incorporating 
temperature in equations to predict AHOD is warranted since 
there may be large temperature differences in the hypolimnia 
of various lakes in summer; Charlton (1980) reported mean 
summer hypolimnetic temperatures ranging from 4.0 to 11.0°C 
in the Laurentian Great Lakes. During the period of 
ice-cover, the temperature variation between lakes is not as 
great as the temperature variation between lake hypolimnia 
during summer; Jackson (1979) reported mean winter 
volume-weighted temperatures ranging from 3.0 to 3.4°C in 
three lakes in Ontario, and the variation in mean winter 
volume-weighted temperatures in the Albertan lakes was 2.2 
to 3.4°C. Although biological processes are influenced by 
temperature, there is little variation in temperature 
between ice-covered lakes, and thus the influence of 
temperature on WODR is likely less than on AHOD. 

The models developed to predict WODR over a broad range 
of lake types need to be tested with an independent data set 
to determine their accuracy. Data from lakes outside the 
original data set, e.g. shallow and unproductive, need to be 


collected to improve the generality of the models. 
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Table 1; Percent soxygendsaturationsing13cAlberta lakes 
during fall and early winter, 1982. Measured oxygen 
Saturation on the last visit before freeze-up (Prior), the 
first visit after freeze-up (Post), and predicted oxygen 
Saturation on date of freeze-up calculated from the 
Y-intercept of the regression of transformed oxygen mass 
against time (Predicted). Numbers in parenthesis refer to 
the number of days before or after freeze-up when the lake 
was visited.f indicates lake did not mix to the bottom at 


fall turnover. 


% OXYGEN SATURATION 


Measured Predicted 

LAKE Pron Post 

Amisk N 64(5) 64(21) 66 
Amisk § 52.6) 54(21) 55 
Baptiste §S bei(5) 6221) 70 
Eden 6201) ee te?) G7 
Halfmoon 48(11) 55 G9!) 64 
Hasse 94(9) Sea) 97 
Hubblest See) 54(9) a5 
Lessard Son sh) 82(14) iM 
Nakamun 60(9) 52.612) 52 
Peanut 66(21) LO) 68 
Saues D2 hela) 72¢40) 7) 
Wnt 47(4) sen) 51 
Wizard 88(9) 67115) 81 
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Table 2: Winter oxygen depletion rates (g O2,:m~?+-d~') during 

the first 112 days past freeze-up, and correlation coefficients, 
pyetorethe mass of oxygen (on lineareandstranstormedpax%= 

K°-** data) vs days past freeze-up: All, rate covering entire 
12-day period; Early, rate based on first three sampling 

dates; Later, rate based on remaining dates; Trans, rate 
calculated on All data transformed. Sample size is indicated 

in parenthesis. t+ indicates lake did not mix to bottom at 


PaliecurnoOvVer. 


WINTER OXYGEN DEPLETION RATES r 

LAKE ALL EARLY LATER TRANS ALL TRANS 
Amisk N 0.554(5) aa aoe OL0877 0.98 0.99 
Amisk S 0.848(5) a = OF 1762 0.98 0.99 
Bace1Stessel. 751s) === a= 0.1180 0.99 0.99 
Eden 0. 33HHw> 30.48.11 3)007286(4)" One 562 0.99 1.00 


Halfmoon OF 6216)" 0537303) 907643843)" (ORS 52 Ome Saas 
Hasse Ore G6) ow 26(3 ) ORNS Ce CRO 52 PCO seiic..0.0 
Hubblest 02475 (7) 002.345(2) N08 520(4) 050846 Oe Sars 


Lessard Om 2468) 0.3421 3a 2225) 00 oTe O.995= 0.32 


Nakamun OF 28010) MO 56743) 30.84 Ges 0.207 73 0.960" 0..98 
Peanut Ore 40S: Gon Ole 12 0S) 06 ona) 08 23) 100) WOO 
Sauer SIRENS) Wi NCD) IAI) Oe Oersis ORGS 9 pale 0 
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Wizard On 5Scs) Zia aie Ore0S2 Weise Whe i)s} 


I 


: iieies. aSniL\ nO: i 
sii s09 unrevo> ated, [iA abe: 28g ae ae Na 


end Lamas gendir cratha jee tis ei aee- toe: — 
soar yanait sag¥ad phinisite? ne besge 976 , tae 
hele vibyiit @i.saia viquae , ee 4 e970 al ao 


$e moagged oc: Abt ton Oe Aist ag: nod tan 4 5 


————s —— 


i 


a 


 ——————$——— 


4 _- BST AR Astras opera \nacite 


242eT 197A 2427 Pars ¥ISaS 324 
EeSG, 88.0 40,0 Sn a \2)a2eq0 
29,0.) Se. et ak an we *)685.8 : 


} 
op. 0) Be IC oe a — (Eka ae 


bom BELO Seen. (aleBS O (eq Bee LIME 
the AO UF Terk r re bow La) 4 ic Vee e sf eae 
ne eel 0, Sept peat) (ERT R “Leos ei 


veh abet easter <pyare lo (ries as | Soke te ibe! 
ge. .aRD apao,o (h1d8e.0 cleslity (Me2TRE 
; a ie gr oa rt i ® {ies 

nO°°A ; ara fae E he Pad 
ao FO’ aT eur } ig}kete oy 
EPTO.U: (7728 101 ENN CE aD RDA 
Bost Jag) fea0.0 (2)e0e00" Verte, 4 

: re 

ot 62,9 9 BBC. 9 1e) 895.0 PEeaeeo 


i 
¥ 
: om | 
G 
o 
«. 
cc 


Cid ah 


43 


Table 32 Correlation coefficients for winter oxygen 
devleciGnecaces=(GrO5emy a hornet ranstoricdn A =.) 
and linear data vs six parameters: mean depth (z), maximum 
depth (Zmax), loss on ignition (LOI), spring and summer total 
phosphorus (TPsp and TPsu), and summer chlorophyll a (Chl a), 
in Albertan lakes. Phosphorus and chlorophyll are expressed 
per unit area of the lake and per unit volume of the 
trophogenic zone. Significant correlations are indicated by: 


ee. 05 ie ee P<0 0 pand +t*(P<02005), 


Parameter Transformed Linear 
z (m) 0.41 0.90%*x 
Zmax (m) Onsi6 0.82*x 
LOI (%) =O 1.24 =05.36 
TPsp (mg-m *) 0.46 0.73%* 
TPsu (mg-m7?) 0.72% 0.64% 
Chi a “(mg-m 7) Ons O33. 
TPsp (mg-m °*) 0.68x 0523 
TPsu (mg-m °) 0.46 On On 


Chl a (mg-m-?) 0.40 -0.09 
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Table 5: Winter oxygen depletion rates (WODR) (g 
O7°m---d7*), ‘mean depth (2) total phosphorus (TP) (mg-m-?), 
and chlorophyll @ (Chl a) (mgq-m-*) from four Sets of 


temperate zone lakes in Canada. 


SPRING SUMMER SUMMER 


LAKE WODR Zz TP TP GHita 


(nm) mgm =)  (mg-m= 2) (mg-m- =) 


ALBERTA ' 

Amisk N 02554 107 902 Zee 61 
Amisk §S 0.848 19.4 422 219 86 
Baptiste S Or TES) Waaie 320 202 64 
Eden OR 3 aes ORIN eal Zul 
Halfmoon On 462 4.8 25), ZAS 147 
Hasse Oro 2 | are Se) 174 26 
Lessard 0.243 Sra 115 oF 42 
Nakamun Os uel) 4.5 181 135 72 
Peanut 0.407 oye 260 Zin oy, 
Sauer ORS eK) 4.2 265 208 28 
Wizard OR oss G2 235 132 64 
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(Table 5 cont'd) 


2s tee ntl rs ae rE ERT eee, OPEN all at ene wees ae ee Maced 
SPRING SUMMER SUMMER 
LAKE WODR Zz TP TP CH1l a 


(m) (mgeim= 2) (mgm 2). (mam =) 


re ee ee 


Four Mile O.28:2 is 124 98 16 
Green Shs aug Gia oF 87 16 
Habiburton ORe2 95 13.6 82 58 12 
Halls On-373 Zilewee We 63 13 
L.Boshkung Oe127 6 80 2) 13 
Maple 0.242 tate 72 79 Te 
Moose Og aS) 16526 wT 13 is 
Oblong 0.186 1 ges 66 Ti 
Pine Ce Ohery, 7.4 1) 85 ilies 
Twelve Mile Onzo6 ies 76 68 16 
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SPRING SUMMER SUMMER 
LAKE WODR Zz TP TP CHl a 
(mm) (mg<m- 7) (mg-m-2) (ing-m- =) 


ee Oe ee ee eee eee 
MANITOBA® ° 


885 0.260 8 184 
882 03270 Ziel 90 
ahs) OS SIE [Bois Sil 
318 Oyaz 20 leat 64 
587 0.260 Pars) 216 
721 0.240 36 26 
675 0.340 mh 55 
200 O30 26 54 
019 03310 3.4 49 
Nora 0.420 4,2 Sy 


"Data on TP and Chl a@ for Nakamun and Halfmoon lakes from 
Riley (1983); all other TP and Chl @ data from Prepas and 
Vickery (1984). 

7Data on WODR from Welch et al. (1976). 

‘Data on Spring TP and summer Chl a@ from Dillon and Rrglex 
(1974). 

*Data on summer TP from Dr. P.J. Dillon (unpublished). 
‘Data on WODR from Schindler (1971). 

‘Data on summer TP from Reid et al. (1975) and G. Linsey 


(Freshwater Institute, Winnipeg, Manitoba, Pers. Comm. ) 
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(Tanke ss. cont'd) 


7Data on summer Chl a from Armstrong and Schindler (1971) 
and (Reid ethal 7UGi975) 

®Data on WODR from Barica and Mathias (1979). 

*Data on summer Chl a@ from Barica et al. (1978) and Mathias 


and Barica (1980). 
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Table 6: Regression analysis of winter oxygen depletion 
rates (WODR) (g O2-m-?-d-') vs four parameters: mean depth 
(z), spring and summer total phosphorus (TPsp and TPsu) 
‘mgem “)) and summer Ghilorophyliea (enl aietma-m 9) for ithe 
lakes in Table 5. Number of lakes used in each regression 
(Al, cCOurelation coefficient of each regression (fF). 
Significant correlations are incicated by *(P<0.05), 


te ( P0005) ,-and **#(P<0.001). 


PARAMETER 1 EQUATION r 

Z a8 WWODR’= Cac l4 + e0700Go4es 0.34% 
TPsp 27) PRODR =e0n 24 = C00 126 iPsp 0.81%**x 
TPsu So WWODR  =20.052 + 0700227 ETPsu O.77*** 
Ciieg 48 “WODR = 0-242 4+ 0.00201 chia 0.44%% 
TPsp,z 27 WODR = -0.001 + 0.00134 TPsp + 0.0119 z 0.89%** 
MPSunzZ 33. WODR = -0.101 + 0.00247 TPsu + 0.0134 z 0.90*** 


Chia 2 438 MWODR = 0.093 + 0.00257 Ch! a + 0.0127 2 O.65"%* 
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Table. 7: Dissolved oxygen (mg-L-') profiles at Six sampling 


Sites on Wizard Lake, March 7, 1983. 
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III. General Discussion and Conclusions 


(1) One-site sampling 

The development of empirical models to predict WODR are 
based on oxygen profiles obtained from a Single site in the 
lake (Welch et al. 1976; Barica and Mathias 1979: Jackson 
and Lasenby 1982). I have shown that one-site sampling is 
not representative of whole-lake oxygen profiles in Albertan 
lakes. Possible causes of this horizontal heterogeneity may 
be differential oxygen depletion (Linsey 1981) or localized 
inputs of oxygen to the lake. Research into the extent, 
magnitude, and causes of this heterogeneity needs to be 
performed to obtain better estimates of whole-lake oxygen 


content under ice-cover. 


(2) Chlorophyll @ under ice 

In 12 of the study lakes the average Chl a level, 
during the period WODR were calculated, was 2.9 mg-m'°*. 
After this period Chl @ levels increased and a concommitant 
rise in DO levels was observed. In the 13th lake, Halfmoon 
bake, "bhere was tan algal bloom shortly after freeze-up (32 
mg-m-* Chl @) and, when the algae died, the WODR for this 
lakeweinereased from 0.373 to 0.448 g O2-m *<d='. ‘Two 
possible explanations for the increased WODR in Halfmoon 
Lake after the algae died are: (1) oxygen may have been 


introduced to the lake via photosynthesis, during the bloom, 
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thereby slowing down the net consumption of oxygen, or (2) 
when the algae died-off, the consumption of oxygen may have 
increased Since the algae were a fresh supply of organic 
material available for decomposition. Since primary 
production was not measured, it is not possible to determine 
the exact effect algae had on WODR in Halfmoon Lake. Thus, 
more research is needed to determine what effect(s) 


7 Ot 


direct and indirect, algae under ice-cover have on WODR. 


(3) Nonlinear WODR 
In this study, WODR were found to be higher at the 

beginning of the winter than at the end. Jackson and Lasenby 
(1982) also found that oxygen was depleted in a nonlinear 
fashion in lakes under ice-cover in Ontario. Transformed 
WODR (WODR') are correlated with estimates of open-water 
productivity, summer TP (TPsu) and Chl a, in the Albertan 
lakes: 

WODR’ = 0.005 = 0700047 TPsu 

WODR" )= 0,062>> 020005Z2eChi ve 


Oxygen will be depleted in a nonlinear manner if any of the 


following conditions are met: (1) the material to be 
oxidized becomes more resistant to breakdown as the winter 
progresses (Hargrave 1971); (2) any oxygen depletion taking 


place in the water column decreases during the winter or; 
(3) the supply of oxygen to the site where oxidation is 
occurring decreases to the point where oxygen supply becomes 


limiting (Mathias and Barica 1980). Since 10 of the 11 study 
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lakes that mixed to the bottom at fall overturn developed 
anoxic conditions at the bottom by late winter, the supply 
of oxygen to the bottom sediments could have influenced 

WODR. The influence which the nature of the material to be 
oxidized and the supply of oxygen to the sediments have in 


DO depletion needs further examination. 


(4) Linear WODR in lakes 

PiSthe WODR in the central Albertan lakes are treated 
as linear for the 112-days past freeze-up, WODR can be 
predicted for these lakes if mean depth is known, an 
observation consistent with the findings of other 
investigators (Welch et al. 1976; Barica and Mathias 1979): 

WODR =90520"> 0208602 
Data from other studies (Schindler 1971: Welch et al. 1976; 
Barica and Mathias 1979) were combined with data from this 
study and a new model to predict WODR was constructed: 
WODR = -0.101 + 0.00247 TPsu + 0.0134 z 

This equation incorporates both morphometry and estimates of 
lake productivity as predictors of WODR, demonstrating the 
interactions of both of these parameters on WODR. This new 
model is the first model which can predict WODR beyond small 
geographic regions with specific lake types. More data need 
to be collected for a wider range of lakes (e.g., shallow 
and unproductive) to improve the generality of WODR models. 
These data can also be used to determine if nonlinear models 


are generally better than linear models to predict WODR. 
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(5) Predicting oxygen profiles under ice 

Existing models to predict oxygen profiles under 
ice-cover (Jackson and Lasenby, 1982) did not accurately 
predict DO levels in the Albertan lakes. The model developed 
for lakes on the Precambrian Shield consistently 
overestimated DO levels while the model developed for lakes 
in limestone basins consistently underestimated DO levels in 
the Albertan lakes. More data on WODR are needed to improve 
on the existing models. Incorporating an estimate of 


productivity will likely improve upon the models. 


(6) Summer AHOD models vs WODR models 

Models developed to predict AHOD (Cornett and Rigler 
1980; Charlton 1980) were able to accurately predict WODR in 
the Albertan lakes. Cornett and Rigler's AHOD models were 
also able to accurately predict WODR over a wide variety of 
lake types. These AHOD models incorporated estimates of both 
productivity and morphometry, illustrating that the 
magnitude of oxygen depletion in lakes is governed by both 


morphometry and productivity. 
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Table 1: Dissolved oxygen (mg-L-'), snow and ice thickness 
(cm) and depth sampled (m) for the Study lakes during each 


VISITE in the fall and winter of 1982-83. 
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Halfmoon Lake 
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Zz 10.41 10.43 8.96 Sie 4.18 2505 2Be24 2.43 
3 NO ee BSCS 2 Crane Seay hc) Sa Sie 2s eA iy al a Ald 
4 iO lcsSe ee ctl 8.47 SOO ohio DiS 20: ee. “Wob2 
5 1039 — 9.19 Gaus 5.40 2609) 2F01Z ees ae eon! 
6 WOR Ses BOOS 7.54 SEs AE ‘eae 1.64 Oso 05 61 
g WOesS Ss leone 4 (Systeie 302 1.56 raed Ole: O7c9" O1720 
8 WORSE. 746 Sesh 2255 O6S 1.29 Oss Cae Rea 
9 O23 Oe Gx ale, e233 aes Oe al 0.88 Oped Oe 20) 
10 OAS lia On OZ Zen6 0.41 0.18 Gass O70Gs 06 
ia EO 2 Ge Ole ets 0.30 OF.03 Ones C2055 Ossie 
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Hasse Lake 


Date Oot 30 Nov 20 Dec 7 Dec 18 waned 
Snow Ts 6 5 5 4 
Ice = ie os 34 45 
Depth 
1 eas Ve 0 11.380 10236 6.09 
Z alee a, ite! 11.64 1 Orie 5.90 
S vias V2gg2 11,46 9.10 Ses 
4 123 Wa 20 10.34 Sag Tago 
5 Teer oat, Zee 3.20 6.02 4,49 
6 ga as Telreedes) oh ANS 6.38 3459 
i) ipiles see) 0256 en Gc 4.75 2693 
8 it29 Te62Z 4,94 1.43 2.200 
9 esi 6.76 Sin Oy C5 fa3A 
Date Feb 6 Feb 16 Mar 2 Mar 14 Mare 7 
Snow 9 9 is 1 
Ice 49 SS: 7) 56 6 1 
Beech 
1 4.82 4.09 3.66 4.26 6s h6 
2 4.46 S56 2.96 2.40 Sas 
3 Aveo Cremns: 1.86 elie 0.96 
4 Shae Zee hee 0.66 0.43 
5 Sj Adee 2.10 1.40 URS C32 
6 Zoe 1.90 O67 O16 0.06 
2 ieee ‘Ake ules 0.54 0.00 O.. 00 
8 ers: Oi 56 Or O:6 COG 0.00 
9 0.00 Ors0i0 O06 0.00 0.00 
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Hubbles Lake 


Date OGne sil Nov 21 Dec 10 Dec 20 Jan 18 
Snow =e 4 2 Gs 2 
Ice re 1) ae S5 oa 
Depth 
1 eos HS 03 6.76 6235 oe ies) 
2 Tad 6.94 6.33 6.33 oa | 
3 TOES Gye | 6.74 6.34 Se 
4 he58 Gmee C#s8 630 Sar 
5 PaOU Cree 62/2 6.85 4.98 
6 7.43 Gags 6.70 6.26 4.93 
a 7.42 6.86 6.66 6.21 4.91 
8 TOO: 6.84 6.49 Gots 4.88 
S) wes 6.51 Gave 6.15 4.64 
10 eke: ere 6.16 6. 03 4.36 
11 Uaeoal G5 heel) 5.66 4.18 
iz leges 6.56 ele as Ses 3.64 
Ve 6.66 6:.5:1 5.60 Bre 3.64 
14 Geos Gas 5.64 5.07 Saal 
es 4.38 6.99 45 4.62 Seeks 
16 05:00 Ces6 5.49 4.19 2690 
Les 02 00 6.26 5.41 4.07 2.78 
18 Or 00 Gy til 5420 SO Lice! 
19 OOO 60) 4,38 $2355 ce 
20 A108, Bs Ol las) Papel 259 
2 1 0.00 5.47 Zee 2.96 0. 00 
Ze 0.00 ‘eae, 1.88 O00 0.00 
23 O00 0.00 0.00 0. 010 0.00 
24 0.00 Grad O 02:00 C00) 0.00 
25 0.00 O00 0.00 0.00 0.00 
26 0. 00 0.00 0.00 0.00 O00 
Oe OOO 0.00 0.00 0.00 0.00 
28 0.00 Gee 00) 0.00 0.00 02.00 
29 0.00 0.00 0.00 0.00 ORAS OLS) 
30 O00 0200 O00 O06 0.00 
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Hubbles lake cont'd 


Date Feb 4 Feb 15 Mar 2 Mar 16 Mat 927 
Snow iD 8 11 9 8 
Bee 54 58 62 65 66 
Depth 
1 ary ope) 3291 3.34 3.43 ee!) 
2 a5 25 Se 2300 210 Seu 
) 4.25 ile! 25.03 2.60 2.690 
a APE) Sie eke Zena Dee 2.68 
=) ho hdl S050 23 2.47 Pe Oe) 
6 4.08 3.49 2a aa0e2 Ua lOrs) 
7 4.07 3.34 2 Oe 202 0.84 
8 S200 SUS 2545 2503 0.40 
2 SAD 2.85 2.04 18 Ona) 
10 Sesiehe) Zee? Zee ear 0.86 
a 2296 22 cel athe 0.44 
ue 2369 1087 eee Cas! 0.46 
be. Peden) laos OnSi6 Ora5S OR 1S 
14 2.20 Se Ol25 0.46 Oe 
15 0G eee Oe 5 Oe exe) O20 
16 legals) Ora 6 0.10 0.138 ORe5 
7 moc Oe eye: OP ONS) O07, Sls his 
18 Ieesi2 0.28 0.00 0.00 0200 
ES) 0.48 O23 0.00 O00 0.00 
20 Oreceuh O.010 0.00 0.00 0.00 
al 0.00 Oi 06 05:00 O00 07.00 
22 O00 0.00 O00 0.00 0.00 
2 O00 0.00 0.00 0.00 0.00 
24 0.00 0.00 0.00 0.00 0.00 
EMS) 0.00 0.00 0.00 0.00 0.00 
26 0.00 O00 0.00 0.00 0.00 
ay 0.00 0.00 0.00 0.00 0.00 
28 0.00 0500 0.00 0.00 0.00 
Zo OF00 0.00 Ou 0.00 0.00 
30 0.00 0.00 0.00 0.00 0.00 
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(Table 1 cont'd) 


Se ae ee eee. 2 te ee ee aie 
Lessard Lake 


Dates Oct 20 Nov 24 pec 11 “pec Za Vale Orca ce 


Snow as 4 5 8 4 4 

Ice ei 25 34 39 48 as 

Depth 

i 1.63 ti. 46 10 S25 Sy ale, Sie el 7226 
2 Wages) eS: 9.95 Oe 7.90 RAE 
3 S204 The ts Sch 8.58 TeOe 6.70 
4 Teks 8.84 Smicus. Sa6e 5.94 Suu 
e Viceia i) gS 4 5u Shee ZG 4.61 
6 de hO 5.96 Bie ehh 4.11 fs3o5 3.19 


Date Feb 2 Feb 10 Feb 23 Mar 14 Mar 27 


Snow 7 7 10 ee 13) 

Ice 56 59 62 65 66 

Depth 

1 Gwo/ 6.46 Si Sere Bae 
2 6.41 5 9 4,95 Sia ae hes TAG! 
S 6.02 Seon 4.15 jbeisies: 1.64 
4 4.85 Sethe Zea ae iis 
3S 4,44 Zoe 2a20 0.99 0.79 
6 1766 1274 1.69 NB 0.28 


Nakamun Lake 


Date Nov 2.) Noy 23° (Dec 1° Dec 11) Dee 22 Dee 29° wan 10 


Snow == 4 3 4 6 4 4 
Ice a 18 24 30 36 46 53 
Depth ba niente oh 
| TS 6252 8.30 (Peuete, 6./0 Smiley 4.59 
Zz Tie Se Sees) WEE Ge 03 SUG 4.69 3268 
3 ise 8.08 eG Seis) 4.69 3.44 eeu 
4 Ve PAs Teele Shoe 4,40 ls oe Lee es 
5 (hvac: 4.73 Zee I} See pete Ceo? Oa 
6 yet Onals 0.96 Oy Behe, ON eis O65 Oe eS) 
¢: Ne Oe ORs 0.42 0.43 O50) Ori Oa22 
8 FOG 0.34 0.24 0.04 G20 Oreriit nS 


Date) dan 22 Feb 2 Feb 10 Feb 23 Mar 14 Mar 27 


Snow 6 6 9 6 6 i 
Ice oi) 62 63 68 pat Us 
Depth 

| 3.44 2206 1.74 Tee Ch WS OS'S 
2 resis leeonill (eye Onay Ores eae 
© 208 Taal nel OFFA 0. 02 0.00 
4 OV 56 1.24 eel 0.49 0.00 Oe OMe: 
5 Om 5.1 0590 G2 56 0.30 On O.00 
6 Oy 4] Sine One OREO? O00 O00 
7 Oats Oe 11 On 0 0500 G00 0200 
8 0.08 O00 0200 0.00 0.00 0.00 
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Peanut Lake 


Date Slelas 255 Nov 19 Dec 8 Dec 18 Uareade. 


Snow OSS 8 rs 6 5 
fce = 12 29 30 39 
Se eee eS A Se 
| 1296 8.84 ie G7 4.71 
2 7.94 8.91 ee) 6665 4.58 
5 7.94 9.03 PeSS 6.03 4.59 
4 7.94 9.04 (ee) 6.44 4.54 
5 esas 8.90 7.34 6u12 4.40 
6 heSe 8.46 62.93 655" Se 22 
7 FP aeie} 7 es 6.14 4.81 2.66 
8 7.84 Wee 4.79 2.44 ie Ss 
9 (nels 6.46 Sete) lice 2 1.34 
10 TSreite) Sieaie) ee Wee.s ORGS 
11 t= 66 Beate iees 1.64 0.34 
12 RAS6 4.50 xo 0 0.61 ORAZ 
13 0500 3.535 0.70 0.44 On02 
14 0 20.0 B50 0.41 02 e0 O00 
Date Feb 2 Feb 10 Feb 23 Mar 14 Mar 27 
Snow 8 11 9 8 a 
Ice 45 48 55 s}e) 60 
Dejeyele 
| Siem OF 230) 2509 Da 
2 Saoe Zarit 1.96 een’! tec 
3 Some) 2.63 1.89 1.339 23 
4 S50s 256 ie eZ 1.09 0.90 
5 2.96 Dee 0.74 OSS 0.44 
6 202 1243 0.86 0.68 Omis 
7 1.04 + 09 O252 O23 0.09 
8 0.86 Oe27 Oess 0.19 Ore 2 
9 O53 Ge26 O45 OR0S 0.06 
10 O32 OraZ O50'5 OG CONG) 
11 Ores O06 0.00 0.00 0.00 
t2 O10 703 O00 O00 0.00 
+3 0.02 0.00 g..00 0.00 0200 
14 0.00 0.00 0.00 0.00 0200 
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(Table 1 cont'd) 


Se a ee eee 


Sauer Lake 


Date Oct 29 Nov 20 Dec 1 Dec 17 Jan 21 
Snow pee 8 6 ' iS 
ice aa 12 19 35 40 
De ee ee en td ee OE ee 
| 8.68 9.46 8.69 On 393 
3 8.61 Sh 8.47 6.86 Sey. 
3 Geos 8.99 sieailts 6.60 Sele 
4 8.49 Gian Tea Ge 6222 2.84 
5 8.49 8.41 7.49 56s Zee 
6 8.49 8.35 tee O10 Bigaed (PAWAS: 
i 8.54 Owes 6.46 4,94 0.89 
8 8.49 a Sipislte: 5be02 0.48 
S) 8.46 1.59 Shas Ses 0247 
LO 8.49 Gmc 4.00 huss Over 
11 Gv62 Geis 2597 Con O2Z3 
te Bre 6 2.93 1.92 0.41 0.12 
1) 8.56 2.46 Gi 7 OS28 0.00 
14 8.50 Pass Oi 42 0.00 02 Oe 
Date Feb 6 Feb 16 Mar 2 Mar 25 
Snow 10 9 14 ie 
Ice 48 a7 48 56 
Depth 
1 2S TKO G20 4.26 
2 ao 0.66 OOS Grace 
3 os 0.44 0.00 Cris 
4 cNesis) ORG O20 GAS 
5 GASC OR10 O20.0 0.00 
6 0.42 0.08 O00 G00 
# Onesie) 0.04 0.00 O00 
8 0.24 O00 0.00 0.00 
9 On19 0. 80 C200 0:00 
10 0.06 0:20:0 0.00 O00 
11 O00 0.00 0.00 0.00 
$2 G00 O'2010 0.00 @<0'0 
ts 0.00 0.00 0.00 OF 00 
14 0.00 0.00 G08 O20:0 
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Twine bake 


Date Nov | Nov 17) Dec 15 Jan 15 Jan 3% Feb 17 Mar & Mar 22 


Snow <= | S y Aes Tal 18 11 
Ice =< 10 37 47 47 Ss) 50 50 
DSR peta ere eae ie SSO cl o> Denia ee SR Ee 
| 9.56 7or46 7.96 Tals Vinod We 1 oe SS 
2 CRE SRs As yehs ees es TEER 1.09 6.03 6.96 
3 Oe Sal MeO 1.67 7.46 fhe 2.6 6297 62600 6.62 
4 SILOS eR vegnacts: Vs iis Toke eee 6.94 Seow ale tuloys) 
5 Rey Men ery Ge hs Ta3e De 3 hep ehe) 6.7/. 6.46 
6 sey palsy”  svegeaal hen THe PHO heal 6.89 6.65 6.45 
7 SRS Uker a teas) | I? Paes 7,09 6.81 Wy opeie ewes! 
8 Shay tomas hi Pas hres 6.96 Binks 6.54 6,25 
9 See tele! Peau eaeS 6.89 Sianis: G0 Su oy, 
10 CeCe WEASIS ee ete) 6.88 G.65 Gry 2s oO 
1 oy Oee ae tom 7.40 7202 6x72 (Seow! 6.19 5.94 
ie OR Oa o6 1.38 (eye SIE 6.49 Shree Galt eee 
i Ciel) ) Vener! fale 6.80 62.81 6.30 Gc UCME os oy 
14 Pe Oe mia 6.97 6 5 6215 6.30 Siete sayeleye} 
15 PLA BS Wi 4 Bile’ 6.88 (Sweet oy! Ein (AS 5.64 5-66) oe 1S 
16 ede 3a G5 Grain 5.647 5. Oe ase) Cones 
17 OFS0seG oS 6.66 Sees Son, ee TAY, Bo) 485 
18 Oris A Texs ie) 6.47 ayrreha| BAG 5) ite ige) Saco 4.76 
19 OP ein acOne2 © Gao 5.70 5.46 Sins) Buco 4.54 
20 Op Olin s .o4 6.25 Sr eh) Sede) Si he 4,76 4.08 
28 Ore Ts oS 5 0 4.45 4.01 4.67 Zee eta 
22 (OOH Oh eae PS 3.94 Beoo Pe IES roe 2. Does oO 
23 OR O07 2027.2 0.95 0.99 leeSul 22S (ec eines) 
24 O00) F 026 0.68 0.47 0.54 0.68 Osu Stee 9 
25 O20 On 207.33 0.34 ORS 0.85 0.04 CRA SO ES: 
26 CO. 00L er 56 0.45 Oh 229) O53 Os Ue. O20 0s 1S 
2] 020057 0.08 OS 0.14 0.04 0.09 C200) O00 
28 G2 007 2.0203 OROF O02 0.04 02.09 O00 On ul 
29 0,00. 50:00 0.00 0.00 0% 0:0 0.00 0200-2000 
30 OU Cs O 200 05 0:0 0200 O00 200 0500 > 0200 
31 0.00, 0.00 0.010 0.00 O00 OF20'0 O00" 462010 
32 eG ae Oe O10) CEC 0200 0500 0.00 02.00) “050 
33 62.007 20'.00 0.00 0.00 0200 0.00 OE OO) eee ee 
34 0200) OF 00 0.00 O00 OF..O0 OF 00 O00 aU Od 
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Table 2: Dissolved oxygen (mg-L-') profiles for the five 


lakes sampled to determine if one-site-sampling is 


representative of the whole lake. 


ee ee ee eee 


Lake Eden, Feb 26, 1983 


Depth 1 2 3 4 5 6 is 


1 Bao SOR, A Ce AES a3 rG aS seo 7s BOuOO 
2 Galo OnO 2 BOT Ou 654 com Scere, o> UG 
3 OPO Gin Se noe On On Garo ea lnmeGe Gt ie asin! G 
4 GeO. <0. OO fo eer oO, 4 aoe Oro OC 
5 DOO met oo 5 Scrat.) 1 e5e08 
6 S629. O02 2404s 7450 
7 202s alin? Simeon 
8 One tl 168 
S OR Size Cree 
10 Cees. “a 
i Oot 02 16 
t2 C203. $003 
a3 0.00 
14 0.00 
15 Coed 
Habimoon Lake, Jan 9, 1983 
SLUTS 
Depth 1 2 3 4 5 6 
1 2 TEP DEO) eS A Oe Boe eA ieee 
Z IiR55> 1e6e le96 02.04 welts 
3 lots Mine S Cel A 6Ow “IO Suir ioe 
4 OMGts Vato Vi. oS. One 4 
5 Oe O) Siete melon wUaIo2 
6 Obvid? wt 60" 20Gy hy 0.64 
7 O29 0750 
8 O.10° 0.04 
Hasse Lake, Feb 6, 1983 
SITE 
SS, ce ed er, I I oS ct eer 
1 GD A160.) 2292 2556 909 Sm0. 00 426 
2 AAG nae eo as 20) S000 
3 i GRE ee ee (aye Sie iis 
4 cn giShay Sic Says | RG IAS 
5 Cpe” Wey oe ee 
6 Pe ee oO ea eT 
q/ boo 
8 O53 
9 0.00 
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(Table 2 cont'd) 


Lessard Lake, Jan 22, 1983 


Srey 
Depth 1 2 3 4 5 6 
] aco Oe0 CH Gln). 2m 2 eo 
Z LU maw ie) Suet) CURRY 16 2aete OG 
& 6770 "6G396 “666e) 6.86 
4 Die 0) S520. wOmor 
= 476) 2.84 
6 o5 ibe) 


Wizard Lake, March 7, 1983 


Depth 1 2 o 4 2 6 

1 Oi ci) Oe Cee oi) Peer 0.co rs nee 
2 OSA SET CAS I, WP ACY. MS ra ots a et prallich ons a18 
3 sry iam Ctrl! S) een UC ee Seale ae Sia 35) 

a ope) eres ks Pes) ere Sh Sis) 

< WSO ee Se ia 2. © Orme ec 

6 Om6 59 2 00ey sO SC 2a 5 

y, OG 9) > 20 | e069 

8 O38" ea 2:0..69 

2) On 2 

10 0.08 

14 O05 
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Table 3: Chlorophyll a (({Chla]) (mg-m-*) in the study lakes 


during the winter of 1982-83, 


ce ee 


Amisk North 


Dec 4 Decree: wanweu 
Depth [Chla] Depth [Chla] Depth [Chla] 
i 13.44 1 2.00 lie 1.48 
2 6.5 ine 1=12 0.46 
3 4.03 3 1.49 (2=15 2On ey 
4 eas, 4 ft nl AS) 20=2607 6" 
Feb 13 Mar 5 Mar 30 
Depth [Chla] Depth [Chla] Depth [Chla] 
3 OF 65 V5 0.54 i=2 C235 
Seals: 0.29 7=9 0.24 6-8 0.24 
O=iSoer Oe OS 20 ee C= SSO nsca 

(Re teal RE lipec ks: i=Ter.0Rs 5 


Amisk South 


Dec 4 Dec 27 Jan 27 
Depth [Chla] Depth [Chla] Depth [Chla] 
| ba 4 1 1.910 1= a OF 62 
2 Sine Z 2.6.0 a ZerOreenG 
3 3.09 3 lesen 14-18 0230 
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(Table 3 cont'd) 
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Lake Eden 
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Depth [Chla] Depth [Chla] Depth [Chla] 
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(Table 3 cont'd) 
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Hasse Lake 
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(Table 3 cont'd) 
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Lessard Lake 
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(Table 3 cont'd) 


a er errr en ot ok ee ee Sea eee EE Sa A et 
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table 2: Mean loss on tenition «(Cl iuobisediment samples, 


and standard deviations (s) of the means, from 12 of the 


study lakes 


Lake LOI S 

Amisk N 39 See 
Amisk § 4] 3.6 
Baptiste §S 29 215 
Halfmoon 30 etl 
Hasse 35 TOem 
Hubbles 38 6.3 
Lessard 63 2.4 
Nakamun 46 39 
Peanut 30 ra 
Sauer 28 4.7 
Twin 30 sts 
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Table 5: Productivity estimates of Hubbles Lake for the 
ice-free season of 1982. Depth of the euphotic zone (EZ), 


chlorophyll a@ concentration ([Chla]) and total phosphorus 


(TP) (mg-m°) an eel es ni ene nena ae Sm 


Date EZ [Chla] TP 
May 14 4 (ise 38 
June 3 3) 2a8 16 
June 24 i) Zins is 


July 15 8.5 3.5 19 


Aug 26 8 CE: 16 
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Table 1: Morphometric characteristics of the study lakes 
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Baptiste South Hubbles 


Surface area=4,429,196 m? 
Volume= 55,856,240 m: 
Maximum depth=25 m 

Mean depth=12.6 m 


Surface area=395,820 m? 
Volume=3,997,871m° 
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(Table 1 cont'd) 
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Lake Eden Wizard Lake 


Surface area=161,460 m? 
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Maximum depth=15 m 
Mean depth=6.95 m 
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Lessard Lake Nakamun Lake 


Surface area=3,207,083 m? Surface area=3,542,528 m? 
Volume=12,492,882 m* Volume=15,768,015 m? 
Maximum depth=6 m Maximum depth=8 m 
Mean depth=3.93 m Mean depth=4.5 m 
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(Table 1 cont'd) 

a ee 
i senessesecsrinrinentunneninenenemsennesenienmnennnee. 
Twin Lake 


Surface area=235,823 m? 
Volume=3,691,504 m? 
Maximum depth=35 m 
Mean depth=15.7 m 
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Figure 1: Bathymetric maps of the study lakes. The sampling 
Sites on each lake are indicated by (m). For lakes where 
more than one site was sampled, the main site is indicated 
by (#1), site number two is indicated by (#2), etc. Only the 


south basin of Baptiste Lake is shown. Contour intervals are 


in Mm. 
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